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§1. Introduction 

Recently, the Calogero- Sutherland mo de# has attracted much interest among many physi- 
cists. The dynamical symmetry of the system enables us to calculate the energy spectrum, 
explicit form of the wave function, and static and dynamical correlation functions. From 
explicit results on these quantities, the system has been identified with a free-particle system 
governed by the fractional exclusion statistics proposed by Haldane.B 

Among the family of the Calogero- Sutherland model, the lattice models such as the 
Haldane-Shastry modeli'S-* and supersymmetric (SUSY) t-J mo deli play important roles 
in the condensed matter physics. These lattice models have the Gutzwiller-Jastrow type 
ground state wavef unction .& In the ground state energy of the t-J model, magnetic and 
charge parts decouple from each other. This has been interpreted as the complete spin 
charge separation.!!' In the study of the dynamics at T = 0, Ha and Haldane identified 
holon, spinon, and antiholon as elementary excitations.^ From these results, we regard the 
t-J model as the fixed-point model realizing the ideal Tomonaga-Luttinger liquid. 

In this paper, we clarify the relation between the t-J model and the fractional exclusion 
statistics from the viewpoint of thermodynamics. For the thermodynamics of t-J model, an 
approach has been proposed by Wang, Liu, and Colemani* and another by Ha and Haldane.!-* 
These approaches, however, rely on empirical assumptions on the state-counting rule. In this 
paper we construct the thermodynamics of the t-J model exactly by using the mapping from 
the continuum Sutherland model in the strong coupling limit 

In the next section, we construct the thermodynamics of the SU(-ft', 1) t-J model by using 
the method of Sutherland and Shastry.0) In this method, we can formulate thermodynamics 
of the SU(i^',l) model for arbitrary value of K. 

In §3-5, we consider the SU(2,1) model. In §3, we show that the t-J model at all temper- 
atures is described in terms of free particles obeying the fractional exclusion statistics.!^ In 
§4, we deal with the low temperature range and reduce the description of the multicompo- 
nent statistics to so-called g-ons. In §5, we present explicit results on the spin and charge 
susceptibilities and the specific heat. These results are interpreted from the viewpoint of the 
g-ons obtained in §4. 

In §6, we extend the low temperature description to the general SU(Jr,l) symmetry. We 
show that the g-on description works well when all the chemical potentials for K species are 
different. On the other hand, if all species share the same value of chemical potential, the 
thermodynamics are described in terms of free parafermions instead of g-ons. 

A shorter account of our results has been presented elsewhere.0) In this paper, we explain 
the calculation in detail and present additional analytical and numerical results. 
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§2. Formulation 

In this section, we construct tliermodynamics of tlie SU(fr,l) t- J model witli 1/ sin^ in- 
teraction. First we explain the relation between the continuous SU(-ft', 1) Sutherland model 
in the strong coupling limit and the SU(i^',l) t-J model. Next we construct thermody- 
namics of the SU(i^",l) t-J model, by using the knowledge of the energy spectrum of the 
SU(-ft',l) Sutherland model. Last we rewrite thermodynamic quantities in a simpler form, 
by introducing a new variable. 

We start from the SU(iir, 1) Sutherland model for particle systems: 

2^9xr L2^sin2[vr(x.-x,)/L]' ^'^ 

where Xi is the coordinate of the i-th particle, L is the linear dimension of the system, and 
A is a coupling parameter. Pij is given in terms of the exchange operator Pij for the internal 
degrees of freedom of particles by 

, —Pa if both i and ?' are fermions 
Pij = { ^ (2.2) 
Pij otherwise. 



The energy spectrum Ex of has been obtainedE^'li^EJ) as follows: 



oo 



Ex= t ^ + ^E E \n'n'\u{nMK') + Eo. (2.3) 



K = 00 k ' = 



Here k runs over integers describing the momemtum. z/(k) is given by z/B(/t) + J2cr=i '^cr('^); 
^'b(k)(= 0, 1, 2 ■ ■ ■) and i'„{k){= 0, 1) are the momentum distribution functions of bosons and 
fermions with species index a (henceforth referred to as spin), respectively. Eq is given by 

^ .'A^iV(iV^-l) 

Let us consider the strong coupling limit A — > oo in order to relate the Sutherland model 
with the J model. In this limit, particles localize with a lattice spacing L/N. Up to 0{X), 
there are two kinds of degrees of freedom; one is the vibration around the lattice points 
and the other is the exchange of particle species between the two lattice points. The former 
corresponds to a free phonon while the latter does to a hopping of an electron or exchange 
of spins. For the Hamiltonian defined by 

n,^, ^ hm f^^^) , (2.5) 



A-^oo \ A 

it was shown!' 00) that 

'Htot = T^t-j + Ti-ph- (2.6) 
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In eq. ( |2.6| ), Hph is Hamiltonian of free phonons: 



2 7 N-1 , 1 ^ 

^ph = ^E?(^-^)(^;^. + 2)' (2.7) 

where J = 2{N/Lf' and (^6^^ is the creation (annihilation ) operator of a phonon. On 
the other hand Tit-j describes the SU(-ft',l) t-Jmodeh 

nt-j = E with A, - ^ sin {"-^^] , (2.8) 

where the subscript i and j represent the lattice points. Fermions with spin a in the Suther- 
land model correspond to hard-core fermions with o", and bosons to holes. In the SU(2,1) 
case, expression (|2.8| ) is equivalent to the following familiar form: 



V. (2.9) 



Here Qo- is the annihilation operator of an electron with spin a; rii and Si are the number and 

J ^ .J 



spin operator, respectively. V = Ylj [I — Cj-^CjiCj^Cjij is the projection operator excluding 
the double occupation. The transfer energy tij and exchange one Jij are given by 2,tij — 
Jij = JD^j^. From now on, we set J = 2 for simplicity. 

We note that thermodynamic quantities of 7it-j are obtained from those of T^tot and ?iph 
by eq. (|2.6|) . In the rest of this section, we construct thermodynamics for Tit-j] first we 
consider thermodynamics for ?-^tot and subsequently subtract the free phonon contribution. 

For the moment, we consider thermodynamics of Titot. Using ( p.3| ) and the relation ( |2.5| ), 
we obtain the expression for the energy spectrum of TitotS'l^^ 

/ rp fP \ ~-2 oo oo 

E,ot = ^lim j = — E E (2.10) 



' K — — oo 



Here we take the thermodynamic limit N oo. In this limit, the energy is rewritten as 

AT POO POO 

Etot = — dk dk'\k- k'\ u{k)u{k'), (2.11) 

oTT J —oo J —oo 

where i/(fc) is to be identified as //(«;) with n = Nk/{27i). In eq. ( |2.10[ ), we see that each 
microscopic eigenstate is uniquely characterized by the momentum distribution functions 
{z/o-(k), z/b(k)}. Hence the entropy Stot per site has the same form as that of the SU(i^',l) 
free particle system: 

1 , 
Stot = I d/c { - V [(1 - z/^) In (1 - v^) + In v„] 

ITX J-oo 

+ (l + z/B)ln(l + z/B) -z/Blnz/B}. (2.12) 
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The density of fermions and bosons are given as 

n^ = — dkiy^{k), = — I (\ki'^{k), (2.13) 

ZTT J-oo ZTT J-oo 

respectively. 

From eqs. (2.12) and ( |2.13| ), we obtain the expression for the thermodynamic 

potentiah Vttot{{v})/N = E^^t {{i^}) /N - Tstot {{i^}) - Ha l^cjn„ {{v}) - ii^n^ {{v})- The 
equihbrium conditions: 

5Vl{{u„, ub})/Su^ = 6Q{{u^, iyB})/SuB = 
yields the equilibrium momentum distribution functions: 

Ua (eik) - fia) = ^ w-rl ^ T ^2.14) 

exp [(e(fc) - ij,^) /T\ + 1 

i/B (e(fc) - /^b) = ^ — ^j7^. -• (2.15) 

exp [(e(fc) - hb) /T\ - 1 

Here we have introduced the one-particle energy e{k) defined by 

eik) = ^^%(K}) = 1 r d^' \k - k'\ u[k'). (2.16) 



N du^ 2 

Equation ( f^.l6D with ( [2. 14] ) and ( |2.15| ) gives the functional equation for e{k). Substituting 
the resultant expression for e{k) into eqs. (|2.11| )- (|2.15|) , we can obtain thermodynamic 
quantities for Titot- However it is tedious to solve eq. ( p.l6| ) directly. In order to simplify 
the calculation, we introduce a new variable p{k) defined as 

p{k) ^ ^ = i 1^ dA;'sgn {k - k') u{k'), (2.17) 

which we call rapidity. In the following paragraphs, we will rewrite the thermodynamic 
quantities in terms of p. 

First we obtain an equation for e as a function of p. Differentiation of eq. ( p.l7|) with k 
gives 

By multiplying both sides of eq. (p.l8|) by p{k) and integrating over k from k to oo, we 
obtain the following relation: 



oo 



eo = / de z/(e 



K 



-Tin (1 - exp [- (e - /ib) /T]) + T J] In (1 + exp [- (e - /..) /T]) , (2.19) 



cr=l 
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where cq = eo(p) = (^r^ — p^)/2 and z/(e) represents v [k (e)). In the derivation of eq. (2.19), 
we have used the boundary condition: 



'pik = ±00) = ±7r, 



(2.20) 



which is obtained from eq. ( p.l7|) . 
For later convenience, we introduce 

X(p) = exp [- (e-/XB)/T] 

and rewrite eq. (2.19) as 

eo (p) = -T In [1 - X (p)] + T In [1 + X (p) exp (C/T)] 



(2.21) 



(2.22) 



(T=l 



with = /io- — /iB. X{p) is the solution satisfying X{p = it) = 0. Equation ( |2.22| ) determines 
X{p) for given p, T, (f^. In the rest of this section, we rewrite the expressions for the 
thermodynamic quantities in terms of p, T, and Co-- Here note that C = /^o- ~ /^b is the 
chemical potential for election with spin a; introduction of an electron is equivalent to 
removing a boson and adding a fermion. 

Now we consider the expression for -Etot- The constraint ub + J2a=i n-a = ^ leads the 
expression for /ie: 



/iB=TlnX(p = 0,r,{C})+ / dp 



TT — p 

u{p) ' 



(2.23) 



where i^(p) is the abbreviation of u {k (p)). We will derive eq. ( p.23| ) in Appendix A. From 
eqs. (|2.21| ) and (|2.23|) , we obtain the expression for the one-particle energy e{p) 

^X{p) 



e{p) = —Tin 



X(0) 



+ 



dp 



n — p 
u{p) ' 



By using ( p.l8|) , (|2.20|) , and ( p.24|) , the internal energy ( p.ll|) is written as 



N 



dpe(p) = — T / dp In 



Using the following mathematical relation: 

' Xjp) 
X(0) 



dp In 

we find a simple form for Etot'- 

27rEtot _ 
N 



1 
T 



X{p) 
X(0) 

p(7r -p) 



dp 



TT — p 

u{p) ■ 



dp- 



dp 



vr^ — p^ 
u{p) 



dp 



u{j>) 
eo(p) 



z/(p) 



(2.24) 



(2.25) 



(2.26) 



(2.27) 



The relation (|2.26|) will be proven in Appendix B. Similarly we can represent density and 
entropy in simple forms. By using ( |2.18| ), the electron density n„ with a spin is written as 



n„ 



1 

2^ 



dp 



z/(p) ■ 



(2.28) 
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The entropy (2.12) per site is rewritten as 



27rstnt = 2tt 



K 



IT - r dplnX(p), (2.29) 

J —IT 



\(7=1 / 

which will be derived in Appendix C. 

So far we have constructed thermodynamics for Titot- In order to obtain thermodynamic 
quantities for Ti-t-j, we should subtract the free phonon contribution from the thermodynamic 
quantities for T^tot- The thermodynamic potential fiph for Hpi^ (eq. (p.7|)) is given by 

27^^]ph/A^ = T r dp In [1 - exp {-eo{p)/T)] + vr^S. (2.30) 

J —n 

The entropy per site Sph and the internal energy i?ph are obtained from the expression (|2.30| ). 
The internal energy for the t-J model is now given by 

27fEt-j ^ 27r(Etot - ^ph) 
N N 

The thermodynamic potential Qt-j is represented by 

2nSl..j _,JE^_ j.^^^ _ I; _ !k) (2.32) 



\ (T = l 



7r3 



= T r dp In ( f , , ) - -. (2.33) 

J-^ \l- exp [-Co (p) /T] ; 3 

The expression (2.33) is a main result in this section and will play crucial roles in section 5 
and 6. It provides a practical procedure to calculate the thermodynamic quantities for the 
SU(2,1) case in section 5. It also gives a starting point of the calculation in section 6, where 
we derive the low temperature description for SU(fr,l) model. 

In summary of this section, the thermodynamic potential for the SU(-ft',l) t-Jmodel ( p.8| ) is 
given by (2.33). Here X(p) is the solution of eq. ( |2.22| ) satisfying the condition X(p = tt) = 0. 

§3. Equivalence to Free Particles Obeying Fractional Exclusion Statistics in 
SU(2,1) Case 

In the previous section, the thermodynamic potential for the 811(^,1) Sutherland model 
were obtained in a simple and closed form. However, it is hard to obtain the physical picture 
of the thermodynamic quantities from eq. (2.33). Hence it is desirable to have an alternative 
description. For the spin model with the 1/ sin^ interaction, Haldanelli-' has already given 
thermodynamics in terms of free spinous obeying semionic statistics. Here, we extend this 
type of description to the SU(2,1) t-J model and construct thermodynamics in terms of free 
spinous and holons. 
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In this section, we first introduce the statistics heuristically, and show that the spinons and 
holons with a particular exclusion statistics reproduces the thermodjTiamics obtained in the 
previous section. Thus the picture of free spinons and holons turns out to be exact. In this 
and following two sections, we represent the electron density by rig = ni + n2, magnetization 
hj m = rii — n2, electron chemical potential by C = (Ci + C2) /2, and magnetic field by 
h^iCi- C2)/2. 

Let us consider a system consisting of three species of free particles: up and down spinons, 
and holons. We require that the energy including the magnetic energy and chemical potential 
is given by 

2vr (^^ - m/i - C^c + C) = 2 / (^TPT + ^iPl + ^hph) , (3.1) 

where = eo(p)/2 — /i is the energy for up spinon and = eoip) /2 + h for down spinon. 
We take eo = eo(p) = (^i"^ ~P^)/2 to make correspondence with the SU(2,1) t- J model. The 
energy of the holon is given by 

eh(p) = C - eo{p)/2. (3.2) 

The distribution functions for spinons and holons are given by Pa{p) and Ph{p), respectively. 

We next introduce the distribution functions for dual particles (holes) by p* for a =|, |, 
and h. Then the entropy is given by 

2^''^ = 9 / dp[(p„ + p;)ln(pa + P*) -Palnpa -p^lnp;] . (3.3) 

The exclusion statisticsi® of each particle is given by the following relation: 

P*a = -Y.9aa'Pa', (3.4) 



where the statistical matrix gaa' has been introduced. Equations ( |3.3| ), (p.4|) give the ther- 
modynamics of free particles obeying the fractional exclusion statistics . 
By taking the functional derivative of the thermodynamic potential: Q/N = E/N — mh — 
— Ts with respect to pa, we obtain the following constitutive equations for the three 
components: 

ea/T = In (1 + Wa) - 9a'a In (l + w~n , (3.5) 



where Wa = Pa/ Pa- From eqs. ( p. 41 ) and (|375|), we can obtain the thermal distribution 
function for each a. In terms of Wa, the thermodynamic potential is written as 

n \ T 



27r 



N 



+ C) = - 2 £ dp E 1^ (l + ^a') ■ (3.6) 
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If we set the following values for the statistical matrix: 





( 1/2 


1/2 


-1/2 \ 


{gaa'} = 


1/2 


1/2 


-1/2 




t 1/2 


1/2 


1/2 ; 



(3.7) 



equations ( |3.5| ) and ( |3.6| ) yield the thermodynamic potential (2.33) with K = 2. Here we 
show it by a direct calculation. 
By using eq. ( ^.51 ) for a =|, we rewrite ( |3.(j| ) as 

27rfi/A^ = -2nC + e| - In (1 + w^) . (3.8) 

With a new quantity X(p) = X{p,T,(, h) defined by the relation: 

exp (eo/T) - 1 



1 + 



expression (|3.8| ) is rewritten as 

2Tin 



N 



T dp In 



X{p)exp[{h + C)/Ty 
Xip) 



(3.9) 



(3.10) 



1-exp [-eo (p) IT] 

Expression ( p.lO|) has the same form as (2.33), if X{p) = X{p). In the following we show 
that X{p) = X{p). 

The constitutive equations ( p. 51) with (|3.71 ) turn into a quadratic equation for Wa] for w^i, 
we obtain the following equation: 

1 _ e^^Vr _ g(2eT+h-C)/T 



Wt + 



1 + e 



(3.11) 



From eqs. ( |3.9| ) and ( p.ll[ ), we can derive the following relation 
eo = -T In 



0. 



1 - X {p)\ +r In (1 + X (p) exp [(C + h) /T]j +T In (l + X (p) exp [(C - h) /T]j . 

(3.12) 

From (|3.9| ), we can see easily that X(p = vr) =0. Thus X{p) is nothing but X{p) and hence 
thermodynamic potential in eq. ( p.lO|) is the same as that in (2.33) for SU(2,1) case. 

Now we find that the t-J model is equivalent to the system of free particles obeying the 
fractional exclusion statistics with the statistical parameters ( p.7|) . This is a main result in 
this section. In the rest of this section, we discuss the meaning of the expression ( |3.7|) . 

Consider an N site system with up spinous, Xj down spinous, and Xh holons. To be 
consistent with the statistical matrix ( |3.7| ), the numbers of available one-particle states for 
1,1, and h should be given by 



D^ = D^ = l + -{N-N^-N^ + N^) 



(3.13a) 



9 



3 



N^-N^-N^). 



(3.13b) 



By the relation dDa/dNis = —g^f^™ we can easily reproduce eq. (p.7|). In the absence of 
holons i.e. A^h = 0, the expression (p. 13a]) reproduces the fractional exclusion statistics of 
spinonsi-*, as it should. We now show that spinous and holons span the full Hilbert space of 
hard core electrons where no double occupation of a site is allowed. With N-^, Ni, A^h fixed, 
the number of many spinon-holon states is given by 



W{N,N^,N^,Ni,) = - 



l + f-l 



[D^ + N^- 1)! {D^ + N^- 1)! (Z^h + iVh - 1)! 



where the factor 



.N-NT--N^-Ni, 



{D^ - ly.N^i (D| - ly.N^i (Dh - i)!A/'h! ■ 

(3.14) 

/2 picks out even numbers of spinon-holon excita- 
tions for even A^, and odd numbers of excited particles for odd A^. We remind that the pair 
creation is an intrinsic property of kink excitations. With eq. ( p.l4|) , the dimension of the 
full Hilbert space is given by 



0<N-^+N^+Ni^<N 



(3.15) 



1/2 \ 
1/2 

1 ; 



(3.16) 



which amounts to 3^. Therefore the spinous and the holons exhaust the Hilbert space for 
the t-J model. 

The statistical matrix given by eq. ( p.7|) is not the unique choice in spanning the complete 
Hilbert space. For example an alternative choice: 

^ 1/2 1/2 

{9aa'} = 1/2 1/2 

^ 

can also span the full Hilbert space. In this case, the holon behaves as a fermion and need 
not be excited in pairs. Obviously the counting of the dimension of the Hilbert space does 
not depend on parameters in a Hamiltonian. The important point is that particles obeying 
the alternative statistics defined by eq. ( |3.1(j| ) are not free in the supersymmetric t-J model. 
Only with the choice eq. ( |3.7| ) the excitations behave as free particles in the model. 

§4. g-on Description for Low Temperature Properties of SU(2,1) Model. 

In the previous section, we constructed the thermodynamics of the SU(2,1) model in terms 
of spinous and holons. In the scheme, spinous and holons are free energetically but interact 
with each other via the statistical interaction. At low temperatures, however, the statistical 
interaction between them plays no crucial roles and, as a result, both spinous and holons 
obey the g-on statistics. The spin-charge separation at zero temperature has been observed 
in ref.(^. In this section, we show that the spin-charge separation still holds up to a certain 
characteristic temperature. 
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Figure (1-a) shows distribution functions of spinous and holons at h = 0. Here = 

1 /2 

' . Both spinous and holons exist only in the region < |p| < vr. However we 
extend the abscissa to the outside of the region, in order to emphasize that the spinon 
distribution function turns into the fermionic distribution function. In Fig. (1-a), we can see 
that spinous and holons distribute in separate regions in the p space. From eqs. ( |3.^j| ) and 
( |3.4| ), the statistical interaction between spinous and holons exist only for the same p. Hence 
spinous and holons decouple from each other in thermodynamics at low temperatures. This 
simple fact indicates the spin-charge separation in this model. At low temperatures, holons 
and spinous can be regarded as g-onslli'lil0), which interact with each other statistically 
only within the same species. 

Figure (1-b) shows the distribution function of up spinous and holons in the presence of 
magnetic field {h > 0). Here p^ is given by (vr^ — Ah)^^^. In this case, few down-spinons are 
excited at low temperatures, since down spinous have a gap in the spectrum. In Fig. (1-b), 
we observe that the spin-charge separation occurs also in the presence of a magnetic field. 
From Fig. (1-a) (1-b), we find that the g-on description depends on the value of h. Hence 
we consider the two cases h = and h ^ separately in the following subsections. 

4-1 Case of h = 

In this case, we define the low temperature region such that the following relation is 
satisfied: 

exp(-C/T)<l. (4.1) 

In the absence of magnetic field, we can set e-^ = = eg, Pi = P[ = Ps, and wi = wi = Wg. 
Eliminating in eq. ( p.5|) , we obtain 

exp(-2es/T) - 1 = (I + w,) exp(-C/T). (4.2) 

Now we define the rapidity region Hg as the region where the following condition is satisfied: 

w,eM-Cm<v, (4.3) 

where 77 is a constant which is much smaller than unity but much larger than exp(— C/T). 
In Hs, the solution of eq. ( ^TTl) is given by 

Ws ~ exp(es/T) for p E H.,. (4.4) 

With the use of eq. ([4.4|), we can rewrite the inequality ( [4.3| ) as 

exp[{e,-0/T]<r] for p G H^. (4.5) 

By definition, the following inequalities are satisfied in the region Hg = [— vr, tt] — Hg comple- 
mentary to Ug'. 

> r/exp(C/T), (4.6) 
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exp{-ejT)<r'exp{-C/T) . (4.7) 
The contribution from the region lis to the thermodynamic potential is neghgible because 

(4.8) 



(4.9) 



0< /_ dpln (1 + w^M < /_ dpln 1 + r/"^ exp (-(/T) < 1. 

Hence the spinon part of the thermodynamic potential becomes 

27r {njN) = -T dp In (l + w;^) - ~^ jj^ ^^^^ + (-^s/^)] • 

Furthermore, we can extend the range of integration in eq. ([4.9|) to [— vr, vr] because 

dpln[l + exp(-es/T)] < / dpln [l + r^"^ exp (-(/T)] < 1. (4.10) 

As a result, we obtain the low temperature expression for the spinon part in the thermody- 
namic potential as 

27r iVljN) ~ -T r dp In [1 + exp {-ejT)] . (4.11) 

From the expression ( [4.11| ), we find that the spinon part behaves as a free fermion one with 
energy eg in the low temperature region. This is because spinous are always excited in pairs 
without participation of holons. 
Similarly, we obtain the low temperature description for the holon part in thermodynamics. 



Wh (1 + Wh) 



We eliminate Ws in eq. ( |3.5| ) and obtain 



exp 



T 



^. (4.12) 

[l_^,exp(-C/T)]^ 

We define the rapidity region lie as the region where Wh < rjexp {(/T) is satisfied. As in the 
case of spinous, main contribution to the holon part in the thermodynamic potential comes 
from the region U^. In He, we obtain 

exp (eh/T) ^ wl/'^ (1 + Wh)^^^ for p eUc- 

Here we introduce Wc = Wc{p) as the positive solution of the following equation: 



(4.13) 



exp (eh/T) = wl/^ (1 + 



a/2 



for p G [— TT, tt] 



(4.14) 



In terms of Wc, the holon part in the thermodynamic potential is rewritten as 



- dpln(l + w^^) - -- dpln(l + w;^). (4.15) 

From eqs. ( ^.13| ) and ( [4.15| ), we find that the holon behaves as a semion with the energy eh 
in the low temperature region. Since the spinon part ( |4.11| ) does not contain the charge 
susceptibility Xc is determined solely by the holon distribution function at low temperatures: 



<9n„ 



Xc 



dC 



h=0 



1 

4:71 



dp 



dph 
dC 



1 

47r 



dp-^ {w, + i/2r\ 



(4.16) 



Equation ( |4.16| ) manifests the spin-charge separation at low temperatures. 
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4.2 Caseofh^O 

In this subsection, we discuss the case where h > and exp {(/T) ^ exp (h/T) ^ 1. The 
holon part takes the same expression as in the previous subsection, while each spinon 
shows a behavior different from each other. Namely, the down spinon does not take part in 
the low temperature thermodynamics because its energy spectrum e|(p) is gapful. The up 
spinon, on the other hand, turns into a free semion with energy e|(p). Eliminating and 
in eq. (|3.5| ), we obtain the following equation for 

1 _ e2VT _ e(2ei-/^~0/T 



^1 + 



1 + e 



-(h+O/T 



0. 



Since wi is the positive solution of eq. ( |4.17|) , we obtain the following relation: 



^(2e^-h-C)/T ^ ^2h/T 



1 



/2> (e2VT_i)/2»i. 



(4.17) 



(4.18) 



From the relation ( [4.18| ), the down spinon part in the thermodynamic potential is estimated 
as 



< 



dp In {l 



w 



C(e 



-2h/T\ 



< 1. 



(4.19) 



Hence it is negligible in the low temperature region. The equation for obtained from eq. 
(P) is 



exp 



2e| 



[w^ + 1- exp {-2h/T)] 



l + iw^ + 1) exp [{h - C) /T] • ^^-^^^ 

As in the previous subsection, we define the rapidity region as the region where the 
following inequalities are satisfied: 

w^exp[{h-0/T]<r,, (4.21) 

w^exp{-2h/T) <r]. (4.22) 

Here we have introduced a constant rj which satisfies the inequalities exp[— — h) /T] <^ 
?7 <^ 1. In H-f, the equation ( |4.2CI| ) becomes 



exp (e|/T) ~ (1 + w-^)^^'^ for p G H|. 
The quantity Ws is the solution of 

exp (cj/T) = wl^"^ (1 + Ws)^^'^ for p G [vr, vr]. 
Then the up spinon contribution to the thermodynamics can be written as 



27rfit 



i.-l£dpi„(i 



W-t 



- J dp\n[l + w~^ 



(4.23) 



(4.24) 



(4.25) 



From eqs. (|4.24| ) and ( [4.25|) , we find that the up spinon behaves as a semion with the energy 
es(p) in the low temperature region defined by eq. ( |4.21| ). 
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§5. Numerical Results 

In this section, we present explicit results for the spin susceptibility Xs, charge susceptibility 
Xc and the specific heat C. First we obtain the integral representation for these quantities, 
by using the results in sections 2 and 3. Then we perform the integral numerically. We 
interpret the explicit results in terms of the g-on picture presented in the last section. We 
will describe temperature dependence of the three quantities, with the electron concentration 
rig and the magnetization m being fixed. Hence C, and h should be regarded as functions of 
ne and m. 

Figure 2 shows the spin susceptibility: 

1 dm 



4 dh 



(5.1) 

/i=0 



as a function of T with rig =1.0, 0.9, 0.8, 0.6, and 0.4 for the electron density. Peak 
structure in Xsi which is a characteristic of antiferromagnetic correlation, can be seen for all 
concentrations. As the density rie decreases from unity, the spin susceptibility Xs decreases 
in the high T region. This simply refiects the decrease of the number of local moments. As 
T — i> 0, the spin susceptibility approaches (2-7?^) ^ irrespective of the electron concentration. 
Even at low temperatures, Xs's for various collapse onto a single curve. Such behavior 
has not been observed in other models such as the Hubbard mode® and the t-J model 
with the nearest neighbor interactionEl'' . Thus we call this behavior "the strong spin-charge 
separation" .0) As temperature increases, thermally excited holon begins to contribute to 
the spin susceptibility and then it begins to depend on the electron concentration. The 
distribution functions for spinous and holons begin to overlap at a temperature where the 
tails of each distribution function have appreciable values at the "Fermi surface" of each 
other. We define the crossover temperature Tx by 

Pc(p = Ps = 7r)=0.01. (5.2) 

We can see from Fig. 2 that Tx serves as a good scale for the crossover phenomena. 
Figure 3 shows temperature dependence of the charge susceptibility: 

for the various electron concentrations. The peak structure in Xc can be seen in the 
low temperature regime for all concentrations n^, = 0.2 ~ 0.9. As T — 0, the charge 
susceptibility approaches 2/[7r^ (1 — rig)]. Xc is enhanced anomalously as the system ap- 
proaches the half- filling, while Xc vanishes just at the half-filling where the charge gap 
is present. This singularity in Xc is a typical behavior near the Mott transition and has 
been observed also in the Hubbard mode© and the nearest neighbor t-J mo deM. Above 
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T ~ 1.2, we can observe the particle-hole symmetry Xc{T,nc = 0.2) ~ Xc{T,nc = 0.8) and 
Xc{T,nc = 0.4) ~ Xc{T,ne = 0.6). This symmetry is expected for localized electrons with 
Xc = "^0(1 ~ nc)/T. Hence we identify T > 1.2 as the high temperature region. 

Now we study to what extent the g-on picture accounts for the low temperature thermo- 
dynamics. The holon is responsible for the charge excitation in the g-on description. In Fig. 
4 the solid curve represents the exact result on Xc which is the same as that in Fig. 3 at 
= 0.8. On the other hand, the dotted curve in Fig. 4 represents the charge susceptibility 
of free holons which is given by Eq. ( |4.1tj| ). The crossover temperature for = 0.8 is 
denoted by x. Below Tx, these two curves are in good agreement with each other. Even 
above Tx, the deviation is small and the free holon picture works well up to T ~ 1.5. Peak 
structure around T ~ 0.06 comes from the van Hove singularity of the free holon band (|3.2| ). 

Figure 5(a) shows temperature dependence of the specific heat: 

1 dEt-j 



a 



t-j 



N dT 



(5.4) 



for electron concentrations Ue =1.0, 0.9, 0.8, 0.6, and 0.4. For all concentrations, peak 
structure can be seen around T = 0.7. Near the half-filling, specific heat changes rapidly 
with temperature in the low temperature regime T < 0.1. Figure 5(b) blows up the specific 
heat in the low temperature region (T = ~ 0.5). The enhancement of 7 = C/T is evident 
for He = 0.9 and 0.8. 

In Fig. 6, we compare the exact specific heat for the J model (solid curve) with the free 
holon contribution (dashed-dotted one), and the free spinon one (dotted one). The sum of 
the specific heat of free holons and free spinous is also shown. The deviation from the exact 
one indicates the contribution from the statistical interaction between spinous and holons. 
The crossover temperature Tx defined in ( ^.2|) is denoted by x. Up to T ~ 0.3, the additive 
contribution from holons and spinous accounts well for the exact result. This confirms the 
spin-charge separation. We can attribute the rapid change below T ~ 0.05 to the free holon, 
which has the small rapidity Pc = it {1 — n^) ~ 0.63 at = 0.8. 

§6. Low Temperature Description for SU(i^, 1) Model 

Here we discuss the low temperature properties of the SU(-ft", 1) long-range t-J model. 
The method in section 3 is not valid for the SU(i^,l) {K > 3) model. Then we follow 
another route; first, we obtain the low temperature expression for the thermodynamics of 
TYtot, and then extract the contribution from the multicomponent J model. We show that 
the internal energy i^tot can be written in the form of the sum of the g-on energy, in the low 
temperature regime defined as the region exp[(„/T) ^ 1 for all a. As in section 4, the low 
temperature description depends on the distribution of chemical potentials for each species. 



15 



In subsection 6.1, we consider the polarized case with Ci > C2 > ■ • ■ > Ck- In 6.2, we consider 
the unpolarized case with Ci = C2 = ■ ■ ■ = C/r = C- 



6.1 Polarized case 

In this case, the low temperature region is defined by the following conditions: 

exp(C/T) > 1 for 1 < a < (6.1) 

and 

exp [(C - C+i) /T] > 1 forl<a<i^-l. (6.2) 

The aim of this subsection is to derive eq. ( |6.32| ). First we consider the rapidity dependence 
of the distribution function z/b(p) and ly^ip) at zero temperature. The distribution functions 
are given by 



00 for \p\ < pc 




TT 



1/2 



for Pc < \p\ 



v„{p) = e{p„- \p\) , with p„ 



cr-l 



TT 



25]Ca + 2(a-l)C 



1/2 



(6.3) 



(6.4) 



Figure 7 shows the rapidity dependence of v{p)~^ , which is a part of the integrand in eq. 
( |2.27| ), for the case of = 2. Near p = Pa (a = c, 1 ~ i^'), //^-(o" 7^ a) can be regarded as a 
constant and hence the rapidity dependence of u'^ is determined mainly by z/^. Hence we 
can replace by 



K 



V 



(7=2 



cr — 1 a— l + z/g- 



K 



1 + wf,i - E 



Gl 



(6.5) 



^ 9o^ = ^/ ~ 1)5 Ga = {a — 1) ^ ioi 2 < a < K, and 



where WF,a = — 1, u^b 
Qc = Gc = K"^. The quantities G^ and Gc give the scaling of the momentum in the g-on 
description. In the particular case of = 2, = 1/2 means that the "Brillouin zone" of 
the holon becomes half of the original size. 

Let us relate the expression (6.5) with the distribution function of g-ons. First we rewrite 
eq. ( p.22| ) in terms of wf,o- and as 

K 



exp(eo/r) = (l + ^i;Bi)n(l 



w 



F.cr 



CT=1 



(6.6) 



We eliminate w-p^^ in eq. ( |6.6| ) with the use of the relation w-p^cj = (1 + tL'B)exp {—(^/T). 
The resultant expression is 



exp(ec/r) = (1 + wb) (1 + w^') [] [1 + (1 + ^b) exp (-C./T)] 



-IjK 



(6.7) 



16 



with ec = {Jlcr Ca — ^q) / K . Here we introduce the rapidity region He where 
i^Bexp {—(k/T) < r] <^ 1 with r] a constant. In He, equation ( |6.7|) reduces to 

exp (ejT) ~ (1 + Wb) (l + ^b')"'^'' = (1 + w^b)'"'^"^ for p G n^. (6.8) 

Then we introduce the g-on quantities pc and Wc as 

Pc = (wc + fl-c)"^ , (6.9) 

exp (ec/T) = (1 + Wc)^^^ (l + w^^Y'^^^ for p G [-vr, vr] . (6.10) 



Apparently, 

(WB+Qcy^^pc (6.11) 

in He. Since both sides of (|6.11| ) are much smaller than unity outside He, the relation ( |6.11|) 
is valid for p G [— vr, vr]. Hence we obtain the following expression for the boson part of the 
internal energy at low temperatures: 

dp — — — ~ / dpeoPc- (6.12) 



Next we consider the term {w-p^a + Qa)^^ in eq. (6.5). With the use of the relations: 
wb = WF,crexp {(a/T) — 1 and wpy = wp.o- exp[(^o- — Co-') /T], we rewrite eq. (|6.6|) as 

exp (eo/T) = (l + w^},) [l - w^^, exp {-(JT)] ' [] { 1 + exp [((.' - (.) /T]} . 

(6.13) 

Here we define the rapidity region as the region where the following relations are satisfied: 

1 < WF,a exp(C/r) for 1 < a < fC, (6.14) 

exp[(C'-C)/T]<«;F,a<exp[(C-C')/^] ioi I < a < a' < K. (6.15) 
In the region Ho-, equation ( |6.13| ) reduces to 

exp (e^/T) ^ w^/^^ (1 + wp.a)^"^" hi 2 < a < K (6.16) 
exp(eo/r) ^ 1 + M;p 1 (6.17) 



with ea = —Co- + Ga (Sct'Ji Co-' ~ ^0 j • Further we define the rapidity distribution function 
for g-ons as p^ = {w^ + QaY^ for 2 < cr < i^, where is given as the solution of 

exp {tjT) = (1 + w^f'^" for p G [-tt, tt] . (6.18) 

As in the boson case (eq. ( |6.12 )), the relation 

dp — - — ~ / dpeoPa (6.19) 



17 



is valid for p G [— vr, vr]. Finally we consider the fermion density: 



27mcr 



dp- 



v{p) ■ 



(6.20) 



As in the case of eq. (6.5), the integrand z/o-/z/ can be replaced by 



V 



Vr, 



a - 1 + z/„ 



K 

E 

(T'=cr+1 



1 



a' - 1 a' - 1 + 



K K + UB 



K 



GaiWa + ga) E ^'a' (^F,<7' + fi-a' ) ^ " (Wr + 5'c) ^ 



(t'=(T+1 



G.P.- E G^,p,, for2<a<ir, 



cr'=(T+l 



cr=2 

Combining eqs. ( PTTFD , ( pTT^ and (|CT|)-(|0^, we obtain 
2 vr (^Etot/iV - Can. + Ci^ 

dpeo 



(6.21) 
(6.22) 



+ 



exp 



K 



Gc / dpccPc + E '^'^ / dpe<^p., (6.23) 



(T=2 



which is valid at low temperatures. The third term in the right-hand side is the contribution 
of g-ons with energy ec, the momentum scaling factor Gc, and the statistical parameter Qc- 
The fourth term can also be regarded as the energy of other g-ons for the spin degrees of 
freedom. 

Similarly, the entropy is expressed by the sum of g-on contributions. First we study the 
boson part of the entropy 

f-TT dp 



/'^ dp 
— [(1 + z/b) In (1 + Ub) - Vb In i^b] 
-n V 



Since v is written as 



K 



^ = w^b' + E [1 + (1 + ^b) exp (-C/T)]-^ 



(6.24) 



(6.25) 



(T=l 



we can replace v by + at low temperatures in the rapidity region He. On the other 
hand, the contribution from the outside of He is negligible (see Appendix D). Hence we can 
rewrite the integrand of the entropy ( 6.24|) as 

(1 + w;B)ln(l + tfB) - welnwB 

1 + Kwb 
(1 + Wc) In (1 -|- Wc) — Wc In Wc 



1 + Kwc 

Gc [(Pc + Pc) In (pc + Pc) - Pc Inpc - p* Inp* 



(6.26) 
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where p* = 1 — QcPc- The relation (6.26) is vahd even outside of 11^ where both w-q and Wc 
are much larger than unity and hence the contribution from the region is negligible. Thus 
we obtain the low temperature expression for Sc as 

/TT 
dp [(Pc + Pc) In (Pc + Pc) - Pc In Pc - Pc In p*] • (6.27) 
-IT 

Similarly, we can rewrite the fermion part of entropy as 



d 

I 

and 



— [(1 — Ui) In (1 — Vi) + Ui In Vi] 

-TT V 

Q rTT 



1 I ^0 

1 - exp ( - — 



(6.28) 



27rs, = -/"—[(!- v„) In (1 - v„) + In v,] (6.29) 
~ a r dp [(p, + pI) In (p, + p:) - p. In p, - pI In p^] for 2 < a < iT, (6.30) 

J — IT 

where p* = 1 — (^o-Pcr- With these preliminaries we can derive the thermodynamic potential 
at low temperatures. From eqs. (6.23), (6.28), and (6.30), we obtain 

2 TT (^Etot/iV - C.n, - Tstot + C?j 

271^3 + T dj9ln (l - e^'o/^) - TG^ J" dp In (l + w'^) -T^Gaf dp In (l + w^^) . 

(6.31) 



(T=2 



The second term in the right-hand side of (6.31) is the same as fiph (eq. ( |2.30| )). Subtracting 
^ph; we obtain the compact expression for the thermodynamic potential ^It-j as 

2n {Qt-j/N + Ci) ^ 71^3 - TG, T dp In (l + w'^) -T^gJ^ dp In (l + w'^) . (6.32) 

The expression (|6.32|) is a main result in this subsection. It reproduces, as it should, the 



result in subsection 4-2 in the SU(2,1) case. 
6.2 Unpolarized case 

In the unpolarized case, we set Co- = C ^-nd ifp.o- = wp for 1 < a < K. Then the low 
temperature region is defined by 

exp (C/T) » 1. (6.33) 

The main purpose in this subsection is to derive the expression (6.40). 
In terms of wb and wp, equation (|2.22|) is rewritten as 

exp (eo/T) = (l + w^') (l + w^')'' . (6.34) 
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The boson contribution is the same as in the previous case; we can replace {wb + Qc) ^ by 
Pc = {wc + gc)~^ for p e [-TT, n]. 

On the other hand, the fermion part is different from the polarized one. With the use of 
the relation wb = wpexp {(/T) — 1, equation ( |6.34| ) is rewritten as 



exp(eo/T)= 1 - Wp^ exp {-(/T) 



-1 



K 



(6.35) 



In the rapidity region Hi where Wp^exp (— C/T) is much smaller than unity, the solution of 
equation ( |6.35|) is given by 



[exp (eo/KT) - 1]'^ for pEU^. 



(6.36) 



With the use of eqs. ( |6.11|) and (|6.36|) , we rewrite i^^tot ( p.27|) . In this case, the integrand in 
( p.27|) is replaced by 



27rE, 



tot 



N 



dp — ^ 



2n^ 
3K 



dpeo 



exp 



1 



1 



KT 



K K + ub 
1 



1 



dpeopc 



(6.37) 



Similarly, the electron density is rewritten as 

/TT 
dpup I V 
-TT 

= 2n -Gc dp {wb + Qc) 

J —TV 

~ 27r - Gr I dppc- 



The entropy is given by 

27rstot - 27rsc 



d_ 
df 



dpT In 



1 — exp 



KT 



(6.38) 



(6.39) 



where Sc is given by eq. (|6.271 ). The second term in the right-hand side in eq. (|6.39| ), 
which comes from the fermion part, is the same as the entropy of free boson with the energy 
eo/K. With the use of eqs. (2.32), (6.37), (6.38), and ( |6.39| ), we obtain the low temperature 
expression for the thermodynamic potential: 



2 TT 



t-j 



N 



TT 



+ n I -I) -TG, dp\n[l + w, 



T I dp In 



3 \K 



+ e 



-{K-l)eo/KT 



(6.40) 



The second term in the right-hand side of eq. (6.40) comes from the boson part and is the 
same as in the polarized case. The third term, which comes from the fermion part, is the 
thermodynamic potential of parafermions with order K, energy eo/K. 
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The low temperature condition ( |6.33|) is satisfied at all temperatures for the SU(ii') 



Haldane-Shastry model, where C oo. Hence the SU(i^') Haldane-Shastry model in the 
absence of external field is thermodynamically equivalent to free parafermions at all temper- 
atures. 

§7. Discussion 

In § 2, we constructed the thermodynamics for SU(-ft',l) t- J model. Here we compare two 
earlier methodJBiB with the method in 5 2. 

One is the method taken by Haldane^-* and Wang et al.B It is based on the empirical rule 
obtained by numerical diagonalization of small systems. When we set /i = 0, the results in 
§4 reproduce those obtained in ref. (|^) as we show in Appendix E. In another limit oo, 
we reproduce the results for the Haldane-Shastry model given in ref. In this approach, 

thermodynamics is given in terms of spinous and holons and hence the low temperature g-on 
description is easily obtained. The trouble with their approach is that it is hard to justify 
microscopically. In principle, we can refer the issue to the Yangian representation theory. In 
the practical sense, however, this approach has not turned out useful in the case of Sl]{K,l) 
{K > 3) case, since the state-counting is too complicated to make an empirical rule for the 
energy spectrum.0) 

Another method has been proposed by Ha and Haldane.i^ They adopted a string hy- 
pothesis for the SU(fC) Haldane-Shastry model. The string hypothesis is applicable to the 
SU(-ft',l) model with a slight modification. However it is not straightforward to obtain the 
low temperature description in terms of free g-ons and parafermions. 

Finally the method0^ we adopted in section 2 provides the easiest way to calculate the 
thermodynamic quantities and to obtain the low temperature description. However it is 
hard to see the thermodynamical equivalence of SU(-R', 1) i- J model with the free fractional 
particles at all temperatures. It has been proven for the special case K = 2 in section 3. 

In the last section, we presented the low temperature description of the thermodynamics 
of t-J model and its multi-component generalization. When all chemical potentials are 
different from each other, the g-on description accounts for the thermodynamics. On the 
other hand, we must invoke free parafermion in order to describe the thermodynamics in the 
unpolarized case. The result in section 6 suggests that a framework unifying the parafermion 
and fractional exclusion statistics should be worth building. 

§8. Summary 

In this paper, we constructed exact thermodynamics of the supersymmetric SU(-ft',l) t-J 



models ( p.8| ) and discussed the relation between the models and the fractional exclusion 
statistics. The thermodynamic potential for the t-J model (|2.8|) is given in a simple form. In 
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the SU(2,1) case, the exact thermodynamic potential is reproduced in terms of free spinons 
and holons obeying the fractional exclusion statistics. At low temperatures, the statistical 
interaction between spinons and holons does not work and the spin-charge separation oc- 
curs. Both spinons and holons obeying the g-on statistics in the temperature region. These 
properties can be seen in the numerical results on the spin and charge susceptibilities and 
specific heat. For the S\J{K,1) { K >3 ) case, the low temperature properties are described 
in terms of g-ons or parafermions, both of which are free particles. 
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Appendix A: Explicit Form for fiB 

In this appendix, we derive the explicit form for /iB- Taking the logarithm of eq. (|2.21| ), 
we obtain 

e = fXB-TlnX. (A-1) 
On the other hand, equation ( |2.16| ) gives 

/■oo 

e{p = 0)=e{k = 0)= / dk'u{k')k' 

Jo 

dpk{p) 

,dk{p') 



TT rp 

dp dp 
Jo op' 

i-TT rp I 

/ dp / dp'—— 
Jo Jo ulp'j 



dp^. (A.2) 

Combining eqs. ( [A - 1| ) and (A-2), the expression for /ib is obtained as eq. ( p.23| ). 

Appendix B: Derivation of Eq. (|2.26|) 

In this appendix, we derive the relation ( |2.26| ). First we prove the following relation: 

d\nX{p) p 
dp Tv{p) 

Differentiating both sides of ( |2.22| ) with respect to p, we obtain 



-p 



^ 1 _i 



- exp(-C./T) +X 1-X 
TudX 
~Y~dp 



dx_ 

dp 
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Tv 



d\YiX 

dp 



(B-2) 



By dividing (B-2) by z/, we obtain the relation Next we integrate partially the left-hand 

side of eq. ( ^ ). 



dp In 



X(0), 



dp (vr — p)' In 

dp (vr — p) 



X{p) 
X(0) 
d\nX{p) 



dp 



^ Tz/ ■ 



From the second to the third lines, we have used the relation (|B-1D . 

Appendix C: Calculation of Stot 

In this appendix, we derive the expression ( |2.29| ). First, we rewrite eq. (2.12) as 



2TTSt 



TT dp 



-n U{p) 



In 



1 + 



nf=i(i-^.) 



K 

+ Y,uM {y^^ - l) + z^B ln(z/B^ + 1) 



cr=l 



From eq. (2.19), we know that 



In 



nli(i-^.) 



T' 



With the use of eq. ( p.21|) , we find that 

In (z/^i - l) = -C/T - InX, + 1) = -InX. 

By using eqs. (|(J-2D and (|U-3|), we obtain 



27rs 



tot 



T 



7^ / dp— 



V 



/ dp 



dp 



(E^=iJ^a + z^B)lnX 



1 
T 



dp 



eo 



dplnX(p)-27r(^C<xn,)/T. 



(T=l 



(B-3) 



(C-1) 



(C-2) 



(C-3) 



(C-4) 
(C-5) 



In going from eqs. (C-4) to (C-5), we have used eq. (|2.28|) . 

Appendix D: Contribution to the Entropy from Outside of He 

In this appendix, we show that the contribution from the outside of the rapidity region He 
to the boson part in entropy eq. (|6.24|) is negligible. 
Since v-q is much smaller than unity outside He, the numerator in eq. (|6.24| ) is rewritten 

as 

[(1 + z/B)ln(l + z/b) - z/bIhz/b] ~ w^^lnw^- (D-1) 
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Except the neighborhood of p = pi = vr, the distribution function z/i is order of unity and 
hence the integrand in eq. ( 6.24 ) is much smaller than unity. In the neighborhood of p = n, 
which is nothing but the region Hi, Ui has the form 



1^1 = (wf,i + 1) ^ ~= 1 - exp (-eo/T) . 

With the use of the relation wb = tfF.i exp {(i/T) — 1, we obtain 

{wBi^y^lnwB ^ (lni(7F,i + Ci/^)exp(-Ci/T) 

^ {Ci/T - In [exp (eo/T) - 1]} exp (-Ci/T) « 1. 

Then we find that the contribution from the region near p = vr is also negligible. 

Appendix E: Reduction to Wang-Liu-Coleman's Result 

Here we show that our result reduces to that of ref. (|^) in the absence of magnetic field. 
Since e^- = = eg with h = 0, we set w-^ = wi = Ws from eq. (|3.5|). In this case, equation 
( p^ ) is written as 



(D-2) 



(D-3) 



ejT = \nws - ^ln(l + Wh^), 
eh/T = In (1 + w;h) + ln(l + w;') - ^ ln(l + w^^). 



The thermodynamic potential eq. ( |3.6|) becomes 

271 {n/N + C) = -T 1" dp In (l + w;') + ^ In (l + w^;' 

From eq. ( |E-1D , we obtain 

^ In (l + w^') + In (l + = In (1 + w,) - ejT. 



dpe^ 



(E-l) 
(E-2) 

(E-3) 
(E.4) 



The substitution O into O yields 

27r {n/N + C) = -T dp In (1 + w. 



T dp In (1 + Ws 



A slight rearrangement of eq. ([E-l|) gives 



Tin 



-1 



Adding both sides of eqs. and ( |E-2D , we obtain 

Tlnwii = ( — Tin {1 + Ws) ■ 

Identifying ( — — TlnWg, and Tlnwh with /x, eg, ec in ref. (J^), we find that eqs. 

( [E-6D , and ( |E-7| ) reproduce eqs. (15) and (16) in ref. (|^). 



(E-5) 



(E.6) 



(E-7) 
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Fig. 1. (a ) p dependence of distribution functions for spinous and holons at zero magnetic field. Zero 
temperature and low temperature results are presented schematically. Here the Fermi surface of holon 
is given by ph = -^/tt^ — 4^. The nonexistent region \p\ > tt is also shown in order to emphasize that 
the spinon distribution function becomes the fermionic distribution function at low temperatures, (b) p 
dependence of distribution functions for up spinous and holons in the presence of magnetic field {h> 0). 

Here Fermi surface of up spinon is given by p| = y/n'^ — Ah. 
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Fig. 2. Temperature dependence of the spin susceptibility Xs for Tie =1.0 (half- filling) , 0.9, 0.8, 0.6, and 
0.4. The crosses on each curve represent the crossover temperature above which the strong spin-charge 
separation breaks down. 



Fig. 3. Temperature dependence of the charge susceptibility Xc for rie = 0.9, 0.8, 0.6, 0.4, and 0.2. The 
crosses on each curve represent the crossover temperature above which the strong spin-charge separation 
breaks down. Peak structure in low temperature regime is due to the divergent density of states at the 
bottom of the free holon band. 



Fig. 4. Temperature dependence of the charge susceptibility for rie = 0.8. The solid curve represents 
the exact result and dotted one does the contribution of free holons. The cross represents the crossover 
temperature Tx . 



Fig. 5. (a) Temperature dependence of the specific heat C for rie = 0.9, 0.8, 0.6, and 0.4. The maximum 
in C around T 0.7 is due to spin excitations. The rapid increase in C for ne=0.9 and 0.8 below T 0.1 
is due to charge excitations, (b) The blow-up of the low temperature behavior of the specific heat. 



Fig. 6. Decomposition of the specific heat into spinon and holon contributions at ric = 0.8. The solid curve 
represents the exact result for the i- J model. The dashed-dottcd one shows the contribution of free holons, 
and the dotted one docs that of fermionic spinous. The sum of the two contributions to the specific heat 
is also shown. Tx is denoted by x . 
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Fig. 7. The schematic profile of i^ip)^^ in the case of -RT = 2. The Fermi surfaces Pc, Pa are shown. Near 
p = Pa, the rapidity dependence of i'{p)~^ is determined mainly by faip) while the behavior of i'{p)~^ 
around Pc is determined by vsip)- 
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